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Abstract 

We use three domain wall simulations from the radiation era to the late time dark energy domination era based on the PRS algorithm 
to calculate the energy-momentum tensor components of domain wall networks in an expanding universe. Unequal time correlators 
in the radiation, matter and cosmological constant epochs are calculated using the scaling regime of each of the simulations. The 
CMB power spectrum of a network of domain walls is determined. The first ever quantitative constraint for the domain wall surface 
tension is obtained using a Markov chain Monte Carlo method; an energy scale of domain walls of 0.93 MeV, which is close but 
below the Zel’dovich bound, is determined. 
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1. Introduction 

All inflationary theories involving Grand Unification Scale 
phase transitions in the early universe predict the existence of 
topological defects HJ. The inflationary field develops a spon¬ 
taneous symmetry breaking phase, where various kinds of topo¬ 
logical defects may form: textures, monopoles, strings, domain 
walls 13 • Although some of these defects leave characteris¬ 
tic patterns in the Cosmic Microwave Background (CMB), no 
conclusive direct observation of defects has been yet found. Of 
these topological defects, cosmic strings have been the most 
studied in the literature, because they were once considered to 
be the primary sources of anisotropy of the CMB Q. The con¬ 
tribution from cosmic strings and textures in the CMB has been 
recently investigated by the Planck probe Q, where stringent 
constraints on these defects were found. 

In the case of domain walls, the situation is different, be¬ 
cause their energy density grows faster than that of radiation 
and matter as the universe is expanding and hence they would 
eventually become the dominant part of the energy of the uni¬ 
verse 151161. For this reason, their symmetry breaking scale is 
believed to be constrained to around 1 MeV. Nevertheless, their 
existence is predicted by various cosmological models which 
have discrete broken symmetries EISl. The purpose of this pa¬ 
per is to carry out a more quantitative derivation of this bound. 

In this paper, we briefly introduce the equations governing 
the evolution of domain walls in the universe. We use the results 
of the simulations to calculate the energy-momentum tensor of 
domain walls for each time considered. We then determine the 


unequal time correlators (UETCs) of the stress-energy tensor 
components at high resolution and precision, relevant for the 
Planck satellite. We use two simulations, covering the entire 
period from recombination to late-time A domination. We di- 
agonalise the UETCs and use the eigenvectors as sources for 
the CMB fluctuations, thus obtaining the domain-wall angu¬ 
lar power spectrum. We combine the results from the three 
epoch considered in order to obtain the final power spectrum. 
Finally, we use a Markov chain Monte-Carlo parameter esti¬ 
mation (COSMOMC) code with the latest CMB likelihoods in 
order to determine the first-ever constraints on the domain wall 
amplitudes. 


2. Domain walls equations of motion 


Domain walls are the simplest cosmological defects, as they 
can simply be described by a single scalar held (p. One starts 
with the Lagrangian describing a discrete broken symmetry, 
that can be written as: 




1 

An 


+ V(<p) 


( 1 ) 


where the held cp is real and the potential V has at least two 
degenerate minima 0. The energy-momentum tensor of the 
walls network can then be expressed in terms of this Lagrangian 
as follows: 
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In a Friedmann-Lemaitre-Robertson-Walker expanding uni¬ 
verse the metric is taken to be: 
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g^y = a^diag(-l, 1,1,1) 


(3) 
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where the 0'^ dimension corresponds to conformal time. With 
respect to this metric, the components of the energy-momentum 
tensor become: 
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where prime denotes a derivative with respect to conformal time 
and the Laplacian is expressed in physical coordinates. By ap¬ 
plying the standard variational technique: 
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for the Lagrange density (jT]) with metric ([^ the equation of 
motion for (p is obtained: 
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where 77 is the conformal time. In the case of a power-law ex¬ 
pansion of the universe, a oc Z and has the value 1 

in the radiation era and 2 in the matter era. 

In order to calculate the stress-energy tensor components one 
has to first solve equation of motion ([^ and then to substitute 
the solution into the corresponding equations Q-(j^. However 
in the comoving coordinates described above, the thickness of 
the walls decreases as and, as Eq. ^ has to be solved nu¬ 
merically on a grid, the wall thickness quickly becomes smaller 
than the grid spacing. This problem can be overcome by modi¬ 
fying Eq. (|^ 0 to: 


5^0 / dlna \ 1 50 2 
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The unmodified equation of motion corresponds Xo a - p - 2. 
However, taking the coefficients to be a = 3 and /3 - 0, one 
solves the problem of wall thinning, as the walls would have 
constant thickness in comoving coordinates (by modifying yS) 
and would also maintain energy-momentum conservation (by 
modifying a as well). The procedure is called the Press-Ryden- 
Spergel (PRS) algorithm, after the names of the authors in 191 ■ 
This equation can now be solved numerically on a grid using a 
finite difference scheme as follows: 
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where dots denote derivatives with respect to conformal time. 
The effects of this modification of the equations of motion have 
been studied in detail by various authors, most recently in IfTOl . 
and no significant effect on the domain wall characteristics has 
been observed. 

These equations use the assumption that the domain walls al¬ 
ways have a small contribution on the overall energy density of 
the universe. This assumption is based on the fact that no direct 
signs of domain walls have been observed. Hence their contri¬ 
bution to the matter perturbations can be treated as a first order 
approximation in perturbation theory. Therefore, at this order, 
their evolution does not significantly affect the expansion of the 
universe and hence we can safely use a power law expansion 
rate for radiation and matter epochs. 

3. Formalism for calculating the power spectrum 

After obtaining the field 0 and its time derivative we use Eqs. 
to calculate the energy-momentum tensor, and then we 
project it onto a 3-dimensional grid as in Ref. IfTTI . This stress- 
energy tensor is decomposed into its scalar, vector and tensor 
parts, and as in the case of the strings, we chose to use two 
scalar components ( 0 oo, ), one vector and one tensor compo¬ 

nents. We have chosen these particular variables, as they appear 
in the matter perturbation equation in the synchronous gauge. 
These get modified by the presence of domain walls as follows: 


l—h‘- 

a 


l-h' 

a 


a ■ ■ 

l-h‘ 

a 


, ^ 4;7rG^„ 

(Pi + Pi) Vi -— 0® 

k 

(14) 

-■ \()nG{a^pI.^ H-©^) 

(15) 

1 \bnG{a^piy + @^) 

(16) 

= \6TTG[a^pZ^ + @^) 

(17) 


where h - h - and 0 ^ satisfies the equation 

20^ - 0, 


qD ^ qD 


d k 

-2 -r 

a 32 


„ © 00 ^ 


^00 


(18) 


Domain walls are active sources, and compared to infla¬ 
tionary fluctuations which are only generated at the surface of 
last scattering, they continuously source the metric perturba¬ 
tions. This makes direct CMB calculations impossible due to 
the huge amount of data generated. Eortunately, one may use 
UETCs of domain walls which contain all the relevant informa¬ 
tion at first order in perturbation theory. These UETCs repre¬ 
sent the 2 -point correlation function of different components of 
the energy-momentum tensor. This quantity has been studied 
extensively and has been used in the computation of the power 
spectrum for cosmic strings nuEiiia. Eor domain walls, in 
this particular gauge, we need three scalar UETCs ((0oo©oo), 
(0'^0‘^), ( 0000 '^)), one vector and one tensor UETCs. In 
fact, we have taken two vector and two tensor components and 
checked that they have the same auto-correlators and that their 
cross-correlators vanish. This is due to statistical isotropy m. 
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As domain walls have a different scaling law compared to 
cosmic strings, the UETCs computed directly are non-scaling 
even in a purely radiation or matter epoch. Therefore, in order 
to be able to use them in a Boltzmann code, we have to fix the 
scaling. In each epoch, the scaling behaviour is achieved by 
considering the quantity: 

C{kTukT2) = —^<©(k,Ti)0(k,T2)> (19) 

\jT\T2 

where 0 corresponds to a generic component of the energy- 
momentum tensor. For the cases considered, these are positive 
definite functions, and hence they can be expressed in terms of 
their eigenvectors and positive eigenvalues as follows: 

C {kTx,kT2) = '^AiVi{kTifvi(kT2) 

i 
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with w; = V^t’, Then one would have to substitute the 
energy-momentum tensor component with 


0 (k, t) ^ ^lTWi (kr) = ^/T ■\fAiVi (kr) 


( 21 ) 


due to the fact that we have divided by Vr in Eq. 19 For the 
vector and tensor cases this can be done straightforwardly, but 
for scalars we have to create a matrix formed by all the scalar 
components and then diagonalise it: 
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The first half of each of the eigenvectors corresponds to the ©oo 
part, while the second one to 0'^. The eigenvalues are com¬ 
mon to both. The resulting eigenvectors are used as sources in 
an Einstein-Boltzmann eigensolver separately for scalars, vec¬ 
tors and tensors and the results are then summed up in order to 
obtain the total angular power spectrum components. 

This procedure works for a simulation accross both the radi¬ 
ation and the matter era if one assumes that the same unequal 
time correlators are valid in both epochs. However this is not 
the case and hence we use the UETCs only in their time range 
of validity. The procedure is described in greater detail in Ref. 
ini. For example, when calculating the contribution from the 
matter era, instead of making the substitution from Eq. (|2T|, we 
perform the following transformation: 


©(k,T) 
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(23) 


and similarly in the case of the radiation era. We go from one 
regime to the other by using a suitable smoothing function to 
avoid discontinuities and numerical problems. The total power 
spectrum is the sum of the power spectra obtained from the two 
eras. Although the procedure is not exact, in ca we show that 
indeed the errors are small. Because domain walls are active 
sources and the fact that they are uncorrelated with the primor¬ 
dial fluctuations, one can calculate their power spectrum sep¬ 
arately from the inflationary one and sum up the results at the 


end. This makes the calculation of the domain wall power spec¬ 
trum easier, as we don’t have to worry about other cosmological 
sources. 

In the cosmological constant era it is not possible to use the 
exponential expansion directly, because the expansion rate does 
not have a simple form as in the radiation and matter eras. How¬ 
ever, it is important to use a late-time simulation, because as 
the growth of the domain walls density is greater than the back¬ 
ground, their power spectrum is going to be dominated by the 
late-time contribution. We propose an approximation to the 
consmological constant epoch by considering the universe to 
be expanding with an effective power law, and we determine 
the power A of an expansion rate a oc f ’ that has the same slope 
today, as tha actual expansion rate. We then generalise Eq. ( p^ 
to include the last simulation as well. We use both the two- and 
three-era calculations to determine the constraints on the do¬ 
main walls, and then show that the change due to the last epoch 
is small. 


4. Simulations 

For our domain wall numerical simulations we use a code 
based on the PRS algorithm 19] with the diagnostic tools in¬ 
troduced in ca. This has been successively parallelised and 
optimised to exploit recent HPC developments; the more recent 
version of the code is described in and a forthcoming pub¬ 
lication will describe further developments. 

In order to solve Eqs. we assume the following 

numerical values of the parameters involved: 


a - 3 

(24) 


(25) 


(26) 


This corresponds to an initial wall thickness of Wo = 10. 



Conformal time 


Figure 1: Scaling behaviour of the UETC in the radiation (blue) and matter 
(green) eras 


The grid spacing is taken to be = 1 and the initial con¬ 
ditions are such that (p is taking random values on the grid 
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between -1 and 1 and ^ - 0 everywhere. We ran 3 three- 
dimensional simulations, with a box size of 1024^ points, one 
in radiation era (d = 1 /2), one in the matter era (A - 2/3) and 
one in in late-tiem A-dominated era. 

The energy-momentum tensor of the domain wall network 
can be evaluated at any time, but we are only interested in the 
scaling regime of the simulation for the calculation of UETCs. 
The regime where the network exhibits such behaviour has been 
investigated in detail in ifTOlfT^ . In Figure[T]we show the peaks 
of the (0oo0oo) unequal time correlators calculated using Eq. 
GE), with each of them centred around the value on the jr-axis. 
This shows how good the scaling is in this regime. 

In Figurej^we have plotted four 2-dimensional slices through 
the domain wall network at different times of the simulation, 
and these show how the network becomes less dense over time. 
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Figure 2: Two-dimensional slices of the domain wall network during the scaling 
regime in the matter era (from left to right and top to bottom - roughly equal 
conformal time steps between the beginning and end of simulation) 


For the A-era, taking into account that the contribution today 
would be the most important, we consider a simulation with A = 
Hoto, where to is the age of the universe. Using the values of 
the cosmological parameters from Ref. El , we take an average 
value of d = 0.95 between various likelihoods. 


In a universe with a oc 
walls can be expressed as: 

P = (l 


f'*, the physical density of domain 
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j without introducing significant errors. In 


A’l _ 


1.93. Using 


Using the units in the simulation, we are obtaining the prod¬ 
uct of the first two terms in Eq. El. The domination of the 
power spectrum by the matter era contribution makes it possi¬ 
ble to safely use A 

the matter era, for the scaling regime we took ^ 

Eq. (27 1 and the fact that the background density of the uni¬ 
verse is given by p = forSE' Hence, the coefficient multiplying 
the power spectrum today is given by: 

|67rGfocr)^ 


(28) 


5. Results 

We have calculated the power spectra for the domain walls 
in the temperature and polarisation channels for the radiation 
and matter era simulations separately and then by combining 
the two simulations together. We have also realised this in the 
three-simulations scenario. 

The results show that the radiation era contribution has a sub¬ 
dominant effect. This was expected, because the growth of the 
density of domain walls over time would mean that their most 
significant contribution is at late times. Indeed, as the matter era 
results completely dominate the power spectra, the errors from 
the procedure of combining the simulations become completely 
negligible. In Figure]^ we have plotted separately the contribu¬ 
tions from the radiation and matter epochs to the power spectra 
in the temperature and polarisation channels showing how the 
matter era dominates on all the scales of interest. Only in the 
scalar TE plot one can see a more important effect of the radi¬ 
ation era. The fact that the power spectra are completely domi¬ 
nated by matter is in agreement to the domination of the results 
by the late-time contributions and also on the necessity to run a 
simulation with a higher expansion rate. 

By adding the late-time simulation, the peak of the curves at 
I = 2 drops by about one third in all four power spectra consid¬ 
ered (Fig. 1^. In that case, the expansion rate of the universe 
is faster than the growth of the domain wall density and hence 
the power is leaking to intermediate scales, keeping an approx¬ 
imately constant integrated power spectrum. As a consequence 
there is only a small change in the CMB constraint on domain 
walls in the two scenarios. 



Figure 4: Comparison between the domain walls power spectra from radiation 
& matter eras and radiation, matter & A eras. In the plots cr/tQ = 1.5 x 10“^ 


6. CMB constraints on domain walls 

We used the COSMOMC code QSl which uses a Markov 
chain - Monte Carlo method to obtain constraints on the al¬ 
lowed contribution of the domain walls to the CMB power 
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Figure 3: Power spectra of domain walls showing the individual results from the radiation (in blue) and matter (in green) eras as well as total power spectra for 
scalars, vectors and tensors in the temperature and polarisation (EE, TE and BB) channels. In the case of the TE polarisation we have plotted the negative pails with 
dashed lines in the same colours as their positive counterparts. In the plots tr/fo = 1.5 X 10“^ 


spectrum. We had to modify the code to accommodate the 
power spectrum from the domain walls. As the domain wall 
matter perturbations are uncorrelated to the primordial fluctua¬ 
tions, their power spectrum can be calculated separately. This 
is very helpful, because although their spectrum would depend 
on the cosmological parameters, the relative change to the in¬ 
flationary spectrum would be small. For cosmic strings it has 
been checked Km 1201 that this variation is less than 20% and 
this is expected to happen for domain walls as well. Domain 
walls are tightly constrained by their TT power spectrum shape 
and hence the parameter variation impact would not be very 
significant. We have used the standard ACDM six-parameter 
model, together with a parameter quantifying the amplitude of 
the spectrum of the walls together with the latest Planck likeli¬ 
hoods. 

We have analysed the radiation and matter scenario, and also 
one involving a late-time cosmological constant epoch. For 
the radiation and matter scenario, we have obtained a con¬ 
straint on the surface density of the domain walls of cr < 
4.22x 10“®kg/m^ (at 95% confidence level), which corresponds 
to an energy scale of formation for domain walls of 0.96 MeV 

m- 

By considering in addition the cosmological constant era, the 
constraints become cr < 3.85 x 10“®kg/m^ and 0.93 MeV. Both 
are in very good agreement with very rough observational con¬ 
straints based just on the anisotropy constraint 6T/T < 10“^, 


which suggest that their energy scale should be less then 1 Mev 
(the original Zel’dovich bound) 0. 

Even though intuitively one may expect the constraint to 
weaken by adding the cosmological constant era (due to the fact 
that there is less power on very large scales), this does not hap¬ 
pen because there is additional power on intermediate scales. 
There are large error for small I and beyond 1 = 10 the inte¬ 
grated power spectra are almost equivalent. 

The values of the ACDM parameters do not shift signifi¬ 
cantly from the standard best fit ACDM Planck values, with¬ 
out domain walls. This is illustrated in Table [T] This is due to 
the fact that the allowed contribution of domain walls is very 
small, because of the different shape of their temperature power 
spectrum. 

Using these values of the energy scale, we have plotted on 
the same graph in logarithmic scale the standard CMB Planck 
power spectrum ifTTll and the domain walls power spectra, nor¬ 
malised at the 95% confidence level for its surface density (Fig. 
1^. The plot shows that indeed the domain walls only contribute 
on large scales as their power spectrum is quickly decaying in 
/-space. 

7. Conclusions 

In this paper we have used high-resolution simulations based 
on the PRS algorithm to evaluate the energy-momentum ten- 
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